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Abstract 

We investigate entanglement between collective operators of two wave- 
packets of finite spectral bandwidth, in two different approximations of 
the Multimode Parametric-Down Conversion (MPDC) process: the pair- 
wise and the one-to-all interaction patterns. For collective operators we 
choose the macroscopic amplitudes of each wave-packet defined by the 
Fourier Transform of their microscopic mode amplitudes. This approach 
intends, to respond to realistic experimental conditions, where measure- 
ments apparatuses may not resolve single microscopic mode amplitudes 
but rather the collective amplitude of the wave-packets. To quantify the 
bipartite macroscopic entanglement we use the logarithmic negativity. We 
relate the time dependent degree of macroscopic entanglement with the 
complexity (number of modes and interaction pattern) and the temper- 
ature of the system. Our results show that the macroscopic entangle- 
ment increases linearly with the number of micro-modes in the case of 
the one-to-all interaction, while in the pairwise interaction it is constant. 
Moreover, in the one-to-all pattern the birth time of entanglement and the 
critical temperature decrease with increasing the number of micro-modes. 
We draw the graphs associated with the two interaction patterns and 
related the degree of collective entanglement with the connectivity and 
the index of each vertex (mode) of the graph. We conclude that quan- 
tum information and computation tasks may be achieved more efficiently 
by manipulating appropriated collective operators in some macroscopic 
systems, then by using their microscopic counterparts. 

PACS number(s) 03.67.Bg, 42.65.Y.j, 03.67.-a 



1 Introduction 



In the last years a great number of experimental and theoretical research has 
been dedicated to find out how to create and extract entanglement from macro- 
scopic interacting systems pQ. Between them are chains of interacting spins 
[21 [3] , harmonic oscillators and harmonic lattices HJ |5J |S] and macroscopic states 
of light [7J [5] . These works intend, to respond to several aspects of the physical 
reality. 

In one hand, in some realistic experimental conditions the measurements ap- 
paratuses may not resolve single constituents composing the systems but rather 
global properties, thermodynamically speaking, they measure macroscopic ob- 
servables rather than microscopic ones. In this case, it is quite pertinent to ask 
wether the quantum correlations existent between the microscopic constituents 
of the systems will survive when the number of these constituents increases? 
On the other hand, thermal states of interacting macroscopic systems, at finite 
temperature, are real classical states and the following question arises natu- 
rally: can interacting systems, initially in a thermal state, evolve to quantum 
entangled states? 

These kind of issues relaunch the old question of the passage from quantum to 
classical behavior with increasing complexity of the systems and with increasing 
temperature [9]. 

Beside this pure fundamental issue, macroscopic entanglement may become 
an important physical resource to perform communication and computation 
tasks with efficiency not achievable classically [TT], with the great advantage 
(comparing with microscopic entanglement) of being created and manipulated 
in simple and accessible macroscopic systems, at room temperature. 

Parametric down-conversion (PDC) processes have been used to visualize 
the quantum to classical transition at the single-mode single-photon level |12j . 
Experimental techniques, developed recently [HI H3], detected quantum entangle- 
ment between bright squeezed vacuum beams (containing about 10 5 photons). 
This suggests that multimode PDC interactions are good candidates to study 
the joint effect of complexity and finite temperature, in macroscopic entangle- 
ment. 

Our perspective in this work will be to identify a small number of collec- 
tive observables (true physical quantities) @], that depend in the size (number 
of microscopic constituents) of the macroscopic systems and such that their 
measurement give information about the degree of macroscopic entanglement. 
We choose for collective observables the macroscopic amplitudes of two wave- 
packets of finite spectral bandwidth, generated in a multimode parametric down 
conversion process, in two different approaches: the pairwise and the one-to-all 
interactions. 

These macroscopic amplitudes are extensive quantities scaling with the num- 
ber of the microscopic modes of the wave-packets, and are defined by the Fourier 
Transform of the amplitudes of the microscopic modes. The corresponding ob- 
servables are macro Bose operators of a two dimensional continuous variable 
Hilbert space. 
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To measure the time dependent degree of entanglement of the bipartite 
macroscopic wave-packets we use the logarithm negativity and analyze its be- 
havior in terms of the size of the wave-packets, and of the temperature of the 
nonlinear medium where the interaction takes place. 

To give a pictorial representation of the physical interactions, we associate 
a graph P3] with each interaction pattern and show that they are fundamen- 
tally different: while the graph of the pairwise interaction is composed of n 
disconnected subgraphs whose vertices have degree 1, the graph of the one- 
to-all interaction is completely connected and its vertices have degree equal to 
the number n of micro-modes in each wave-packet. This structural difference 
between the two graphs will be reflected in the behavior of the degree of macro- 
scopic entanglement. We show that more complex is the system stronger is 
the macroscopic entanglement and higher is the critical temperature where it 
disappears. Another interesting conclusion is that, at finite temperatures, the 
minimum interaction time needed for the wave-packets to be entangled decreases 
with increasing their complexity. 

Our work extends previous studies with bright beams in three aspects: (a) 
our concept of a macroscopic system resides in the number of interacting mi- 
croscopic modes (degrees of freedom), rather than on the amount of photons 
produced during the interaction, and is valid for any finite number of micro- 
modes in each-wave packet, (b) entanglement between macroscopic systems is 
detected through collective operators that retain the global properties of the 
bipartite Gaussian, continuous variables macro-modes rather than computing 
multimode entanglement, (c) the effect of the temperature of the nonlinear 
medium on entanglement is taken into account, previous results only consider 
the initial vacuum state. 

The plan of the paper is the following: In Sec. 2 we introduce the Hamiltonean 
of the nonlinear interacting wave-packets in two different approaches and derive 
their dynamical behavior. In Section 3 we define the time-dependent collective 
amplitudes in terms of the microscopic ones. In Section 4 we compute the 
macroscopic degree of entanglement and the birth time of entanglement in terms 
of the number of micro-modes and of the temperature. Finally in Section 5 we 
present our conclusions. 

2 System dynamics 

In the broadband down-conversion process |15j the pump wave 0, the signal 1 
and the idler 2, behave like short wave-packets of bandwidths Awo, Au)\ and 
Aw2, respectively, interacting in a nonlinear medium. 

In a discrete mode description of the electromagnetic field (e.m.) [16] each 
wave-packet j, is composed by n = 2m + 1 discrete monochromatic modes of 
frequency uij^ = Qj + kS ( j = 0, 1, 2 ; k = — m, — m + I, m — 1, m) where Qj 
is the central frequency of wave-packet j and 8 — is the frequency spacing 
between two neighbor modes. 

We assume that the energy conservation condition is obeyed by the central 
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frequencies of each wave-packet, i.e., (I>o = Q\ + Q2 and we also assume that 
the spectral bandwidth of wave j, Aluj = U)j t m — ujj.- m = 2m6, is small enough 
to prevent overlapping of the three waves and such that each mode interacts 
nonlinearly with all the others, except for those belonging to the same wave- 
packet. The quantum version of the classical Hamiltonean jTTJ [TB] describing 
this nonlinear interaction is 

2 m tci m m 

Hnl = Y E nuj J- kh ],k^],k +ihg Y E E ( a °,?> a lfc4,i ~ h - c -) W 

j=0 k=—m p= — m k=—m l=—m 

where the coupling constant g is proportional to the second order susceptibility 
of the medium. 

When the incident pump wave, is intense, the ensemble of modes (0, k) com- 
posing the pump wave-packet can be treatead classically as a coherent unde- 
pleted field of complex amplitude a — \a \e~ l2v and an arbitrary pump phase 
if. Within this approximation we obtain the following one-to-all Multimode 
Parametric Down Conversion (MPDC) Hamiltonean, 

2 m rn m 

^ = E E to>j,k*i,k*i,k+™* E E ( e " i(<Soi+v) < fc 4,* - h.c) (2) 

j — 1 fc— — m fe=— m l=—m 

where each mode (l,fc) of the signal wave interacts with every mode (2,1) of 
the idler wave and vice-versa. In this approximation the coupling parameter 
w = g\ao\ is proportional to the amplitude of the pump and t is the interaction 
time that is proportional to the length of the propagation path of the e.m. field 
in the nonlinear medium. 

When Acj << Awi, Aw 2 , MPDC Hamiltonean can be further simpli- 
fied, since the coupling between the energy non conserving modes can be ne- 
glected and each mode (1, k) of the signal interacts uniquely with mode (2, —k) 
of the idler, such that u>i k + w 2,-fe = ^o- Within this new approximation the 
Hamiltonian Hi reduces to the pairwise MPDC Hamiltonian [19l [20] 

m 

H 2 = (3) 

k=—nt 

and 

H k = Tiu)i^a\ k ai^ + ^2-ka\ _ k a2- k + iw%(e~ l ( u ° t+v ^ a\ tk a\_ k - h.c.) (4) 

is the Hamiltonian of each pair of field modes interacting under the PDC ap- 
proximation [21 j - 

A better insight about the type of interaction between the microscopic modes 
of the two wave-packets can be achieved when a graph is associated to the 
interaction Hamiltonians H\ and H2, as we do in Fig.l. In panel (a) graph 
G\ = (V, Ei) represents the interaction hamiltonean Hi and in panel (b) graph 
G 2 = (V, E2) is associated with Hamiltonian H.2- In both graphs each vertex 
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(j; k) V = {(j, k) : j = 1, 2 and k = — m, — m + 1, m — 1, m} represents a 
microscopic mode of frequency Uj^ and each edge {(1, fc), (2, Z)} g E represents 
the interaction between mode (1, k) of wave 1 with mode (2, 1) of wave 2. The 
coupling strength w is the same for all edges. 

The degree of vertex (j, k) is the number Afj^k of its neighbors [Hj. In graph 
G\, Mj,k = n, increases with the number n of modes in each wave-packet and, 
even though micro-modes of the same wave-packet do not interact directly, a 
common micro-mode in the other wave-packet is an interaction mediator be- 
tween them, therefore, there is at least a path connecting any pair of vertices 
{(j, k), (i, I)}, in graph G±, i.e., G\ is a connected graph. In graph G 2 , Afj t k = 1 
for any vertex and Gi is composed by n similar, disconnected subgraphs. These 
differences between the degree of the vertices and the connectivity of the two 
graphs must be reflected in the behavior of the degree of entanglement of the 
two wave-packets as we will show. 



3 Collective Operators 

When measuring wave-packets composed by several micro-modes, the measur- 
ing apparatus may not resolve and detect the individual mode amplitudes but 
rather the macroscopic amplitude of each wave-packet. In this case, choosing the 
collective operators as the physical quantities actually detected by the measur- 
ing apparatus, will enable us to relate the degree of macroscopic entanglement 
between wave-packets with the number of their microscopic constituents. 

The frequency dependent collective operators associated with the macro- 
scopic amplitudes of wave-packets 1 and 2, at the initial time t = 0, are nat- 
urally defined by the Fourier Transform of the (n = 2m + 1) modes of each 
wave-packet 

^ tci rn 

AlO = , ^2 ® 10k ' ^20 = / o , , ^2 «20,fc ( 5 ) 

V2m + 1 , ' V2ro + 1 , 

K=—m k — — rn 

the subscript zero refers to initial time. These collective amplitudes obey the 
commutation relations [Aio,Aj ] = Sij and i,j — 1,2. 

The time dependent collective amplitude operators are simply given by 
the dynamical evolution of the initial amplitudes defined above, i.e., Aj(t) = 

AjoUt j = 1,2 where Ut is the unitary time propagator of the composite 
system, that preserves the canonical commutation relations at any time t. 

Let us stress that the collective amplitude operators Aj(i) are true annihi- 
lation operators of the macro-modes (wave-packets) 1 and 2. The macro- 
scopic quantum number operators associated with these macro-modes are nat- 
urally given by Nj — A^Aj and their eigenvectors are collective Fock states 
3 = 1,2. 

Let 

Qj = ^(4+^]) 5 Pj = -^(A i -A\) (6) 
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be the related quadratures phase collective operators of macro-mode j. 

Grouping together the operators Q- } and Pj in a vector of operators X = 
(Qi, Pi, Q2, P2) T , then the canonical commutation relations for the X- s at 
any time t are encoded in the symplectic form [jQ,.Xj] = iOij, where fijj are 
the elements of the symplectic matrix 



e 

k=l 



1 

-1 



(7) 



The initial and time-dependent quadratures phase collective operators are then 
related by the symplectic structure of the dynamical transformation. 

The time dependent amplitudes of the micro-modes, in the pairwise interac- 
tion, are given by eqs. (22) and (23) of Appendix A. Then, the corresponding 
time dependent collective amplitudes have the simple form 



cosh wt \ -> 

Mt)= ^—^ E 

k——7n 



% /2^T 



sinh wt 



V2 



-+ , „ sinhwi v— * 

A(t) = -== E 

V2m + 1 , ' 

k——rr 



iU2,-kt p. 



m 



cosh wt 



i S 



-iuji-kt *t 



k=—m 



V2r 



1 



E 



e a 

„iU2,-kt 



20. 



(8) 



(9) 



In the case of the one-to-all interaction, the time dependent collective ampli 
tudes are 

2m+l 

A(t)=At(t) = E rnj{t)a j0 

2m+l 
3=1 



4m+2 

E n i(*$o 

4m+2 

E ^-(*) s io 

j=2m+2 



(10) 



(11) 



where we used the time dependent amplitudes of the micro-modes, eqs. (30), ( |3l[ ) 
derived in the Appendix A. Note that the Bose operators Ojo,fc were renamed 
for commodity of calculations. The expressions of the coefficients m,j(t), %(<), 
tj(t) and Uj(t), are given by equations (35) and (36) of Appendix B. 



4 Degree of macroscopic entanglement 

In this work, we are interested in the amount of entanglement we can extract 
from the system if only the averaging collective observables Q\ t 2 and Pi,2 were 
measured and manipulated. In particular, we compute the degree of macroscopic 
entanglement in terms of: (a) the number n of modes constituting each wave- 
packet, (b) the interaction time t, (c) the strength of the coupling parameter w 
and (d) the equilibrium temperature T. 

Two physical situations, described by two different initial quantum states, 
have special interest when considering macroscopic correlations and entangle- 
ment in multimode parametric down conversion processes. They are the spon- 
taneous emission and the pair production at temperature T. 
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In the MPDC spontaneous emission each mode (j, k) composing the two 
wave-packets is initially in the vacuum state |0)jo,fc and the initial state of the 
composite system is the pure and separable state given by the tensor product 

iv>>o=®j 2 =1 ®r= 7 jow 

When the nonlinear medium where the waves propagate is in thermal equi- 
librium at temperature T, the initial state of the wave-packets is a separa- 
ble mixed state given by the density operator p = ®j =1 ®^L_ m pjo.fc where 
Ajo.fc = Y^ k =o e_,9j ' fc " j0 ' fc l n jo,fe)(^jo,fe| is the thermal field density operator of 
mode (j,fc), with /^^ = ftwj,fc/ '(ksT) and ks is the Boltzmann constant. The 
initial mean number of photons in mode (j, k) is fijo,k = i/{e^- k — 1). At room 
temperature and in the optical part of the electromagnetic spectrum << 1> 
however, in the microwave part of the spectrum fij^ >> 1 and we cannot ignore 
the presence of thermal photons in the nonlinear medium. Both initial states 
we are considering are Gaussian states. 

Since the Hamiltonians H\ and H 2 are bilinear in the field macro-modes, 
the overall output state is also Gaussian and the separability criterion [551 US] is 
completely characterized by the first and the second statistical moments of the 
collective quadrature phase operators, which will be denoted, respectively, by 
the vector of first moments X = ((Qi), (Pi), (Q2), (P2)) and by the covariance 
matrix (CM) X of elements 



E tf = - (X t X, + - (X t ) (X 5 ) (12) 



The first moments can be adjusted to zero by local unitary operations which 
leave invariant entropy and entanglement. 

The time dependent covariance matrix E(t) for the collective observables 
can be expressed in terms of the three (2 x 2) block matrices a(t), j3(t) and j(t) 

m - ( 3% tl ) (13) 



Yd) m 

The diagonal blocks a(t), j3(t) are the local CM of wave-packets 1 and 2, re- 
spectively. The off-diagonal block, j(t), encode the intermodal correlations 
(quantum and classical) between the two wave-packets. 

The entries Ey, in the one-to-all interaction, computed in the Appendix B, 
(eqs.(|37| - ( |45[ )), depend on the number n of micro-modes in each wave-packet. 
Making n = 1 in these equations we obtain the entries of the macroscopic CM 



for the pairwise interaction, given by eqs.(46)-(54), which are independent of 
the number n of micro-modes. This result is not surprising, since we have 
seen that the graph G2 associated to this interaction is simply a collection of 
n disconnected subgraps, each of them representing the well known two-mode 
PDC process. 

Hamiltonians that are quadratic in the bosonic operators, possess universal 
quantum invariants, i.e., certain combinations of variances which are conserved 
in time independently of the concrete form of coefficients of the Hamiltonian 
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pHl |2"5] . These invariants exist due to the symplectic structure of the transfor- 
mation relating initial and time-dependent values of the quadrature components 
operators. Concerning the statistical moments, these universal invariants are 

m 

X x = det 53 (i) ; X 2 = det a(t) + det f3(t) + 2det 7 (t) (14) 
Defining the quantities 

S(t)=S + i(det7(i)-|det 7 (t)|) ; S =l 1 -]l 2 + ~ (15) 
2 4 lb 

the entanglement criterion can be written in a simple form 

S(t) < (16) 

For both vacuum and thermal initial states and for both interaction patterns, 
we prove in the Appendix B, that 

S(t) = So -\ det j{t)\ (17) 

with 

S = N 10 N 2Q (N W + 1){N 20 + I) (18) 

where Nio and N 2 o are the initial average number of photons in each wave- 



packet, (see eq.(55|). For the initial vacuum state, iVi = iV 2 o = and So = 0. 



4.1 Birth time of entanglement 



Using the entanglement criterion (16), we conclude that the two wave-packets 



are entangled for interaction times ti obeying 

\ det j(U)\ >S (19) 

where |det7(ij)| is an increasing function of t. The instant tg, such that 
I det 7(£f)| = So is the minimum interaction time needed before macroscopic 
entanglement appears, for obvious reasons it will be named the birth time of 
entanglement (BTE). 

For the vacuum initial state and for both pairwise and one-to-all interactions, 
S(t) > for t > 0, therefore the two wave-packets are entangled since the very 
beginning of the interaction. On the contrary, for an initial thermal state, the 



condition (19 1 is attained only after a finite time of interaction, ti = tg, as 
passed. 



For the pairwise interaction the birth time of entanglement 

I 



i I 9 1/2 1 

t £ = — ArcSinh <^ 2^ ^ K (20) 

2w \ N w +N 2a + lj V ' 

is independent of the number n of micro-modes. It depends on the temperature 
T, through the initial average macroscopic number of photons. The dependence 
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of de dimensionless BTE, T£ = wts, given by eq.(20), over the temperature 
is displayed in Fig. 2, for the dimensionless frequencies u>i = 200, Qq, = 400 
and dimensionless spectral bandwidths Au>\ = AH12 = 0.02. We observe that 
BTE increases with the temperature although this growth slows down at higher 
temperatures. 

A similar behavior is found in the case of the one-to-all interaction, see Fig. 3, 
where the BTE, t£, is displayed as a function of the temperature T for different 
number of micro-modes. We observe that for a given n the BTE increases with 
the temperature. We can understand the behavior of the BTE with changing 
temperature in the following way. The MPDC interaction produces correlated 
pair of photons in waves 1 and 2 and entanglement between the wave-packets 
starts only when the system attains the break-even, this is, when the average 
energy of the generated pair exceeds the average energy of the initial thermal 
state. 

The effect of the number of micro-modes in the birth time of entanglement 
is displayed in Fig. 4 for different temperatures. We observe that, for a given 
temperature T, the birth time of entanglement decreases with increasing number 
of modes in each wave-packet. 



4.2 Logarithmic negativity 

To quantify entanglement between the macroscopic wave-packets we use the 
logarithmic negativity j^lj given by 

£/V = max[0, — In v] (21) 

where v is the symplectic eigenvalue of the partially transposed Covariance 
Matrix S. 

Since the CM in the pairwise interaction is independent of the number n of 
micro-modes and since it coincides with the CM of the one-to-all interaction 
for n = 1, we conclude that the macroscopic degree of entanglement £jv, in 
the pairwise interaction, is n independent and its value is the same as if there 
were a single pair of modes in each wave-packet, provided that their frequencies 
are equal to the central frequencies of the two wave-packets and their average 
number of photons are equal to the average number of macroscopic photons in 
each wave-packet. 

On the contrary, in the one-to-all interaction, the complexity of the in- 
teraction pattern is revealed on the macroscopic degree of entanglement. We 
compute £/v as a function of the number n of micro-modes in Fig. 5, for the 
initial vacuum state (triangles) and for a thermal state at equilibrium tempera- 
ture T = 30°K (circles) and for the dimensionless interaction time Tj = 0.3324. 
The same dimensionless frequencies as in Fig. 2 are assumed. The solid lines are 
parallel straight lines. We conclude that for, a given temperature, the degree 
of entanglement increases linearly with the number of micro-modes n, in each 
wave-packet provided that the time of interaction is bigger than the BTE cor- 
responding to that temperature. The value of the slope of the straight lines is 
an increasing function of the coupling parameter w. 



9 



When the initial temperature of the nonlinear medium rises, there are ini- 
tially more uncorrelated thermal photons in the wave-packets. We expect intu- 
itively that, for a given number of micro-modes and for a given interaction time, 
any entanglement should vanish at large enough temperatures. To exemplify 
this behavior we compute the degree of entanglement £j^, for the it one-to-all 
interaction in terms of the temperature T, assuming r, = 0.6978 and for: n = 1 
panel (a), and n = 5 panel (b). The same dimensionless frequencies as in Fig. 2 
are assumed. We observe that entanglement decreases with growing tempera- 
ture until it vanishes at some critical temperature T c — A.38°K for n = 1 and 
T c = 1200° K for n = 5. As mentioned above, the CM for the pairwise pattern 
coincides with the CM of the one-to-all pattern for n = 1 and we conclude 
that, for a given interaction time, the critical temperature T c is independent 
of the number of modes in the pairwise case and increases with n for the one- 
to-all case, contradicting our initial intuition that, more macroscopic are the 
wave-packets more classical should be their collective correlations. 

To connect the physical birth time of entanglement tg and the degree of 
entanglement £j*j, with the coupling parameter w, we use the dimensionless 
time t = wt. Let us assume that for a given temperature T the dimensionless 
birth time of entanglement is Tg, then when the coupling constant w increases 
the real birth time of entanglement tg, decreases, meaning that stronger is the 
nonlinear interaction faster the break-even of energy is attained. In a similar 
way we show that the real interaction time U needed to attain a given value of 
£jV, decreases with increasing w. 

4.3 Interpretation 

At a first glance the growth of the macroscopic degree of entanglement with 
n and the decreasing of the BTE with n are quite surprising behaviors. We 
would expect to see emerging a classical behavior with increasing number of 
constituents of the wave-packets, this is, with their passage from microscopic to 
macroscopic systems. Instead, the collective quantum behavior is insensitive to 
the number of micro-modes in the pairwise interaction and is enhanced in the 
one-to-all interaction. 

This behavior can be explained in the following way. When two micro- 
modes are added to each wave-packet (symmetric towards the central frequency 
uij) that contains initially (n — 2) modes, the interaction energy of the system 
increases roughly as 2nhw (2n terms must be added to the interaction Hamil- 
tonean) in the one-to-all interaction, and increases with 2hw in the pairwise 
interaction. However, the energy contained in the modes themselves (noninter- 
acting Hamiltonean) increases with 2h{ui + ui 2 ) in both interaction patterns. 
The ratio between the interaction energy and, the non-interacting one, increases 
with n in the pattern G\ and is constant in the pattern G 2 - An immediate re- 
sult of this energy balance is that the break-even of energy (BTE) is attained 
earlier and the correlated pair production, that is responsible by the degree of 
entanglement, is enhanced by a factor of n in G\ and keeps constant in the 
pattern G 2 - 
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This is expressed in a simple way by the degree of each vertex in the graphs 
G\ and G2. As a matter of fact, at a given temperature, the degree of entan- 
glement for the one-to-all interaction is a linear function of the degree n of any 
vertex of graph G\ . In the pairwise interaction the degree of any vertex of G2 
is 1, independent of the number of vertices in the pattern, such as the degree of 
entanglement is insensitive to the number n of micro-modes. 

In conclusion, the number of constituents of the macroscopic systems is not 
per se the determinant factor when computing their degree of entanglement, 
the determinant ingredient is no doubt the type of interaction between their 
constituents. 

5 Summary and concluding remarks 

On the fundamental side our work demonstrates that, in principle, purely quan- 
tum correlations can be detected, at finite temperatures, by measuring a few 
number of appropriated collective obscrvablcs of macroscopic systems. 

To illustrate this behavior we studied the macroscopic quantum entangle- 
ment of the collective amplitudes of two wave-packets of finite spectral band- 
width, generated in a parametric down conversion process in two different ap- 
proaches: the pairwise and the one-to-all interactions. 

We have shown that, in the pairwise pattern, the degree of collective entan- 
glement does not depend on the number of micro-modes in the wave-packets and 
that, in the one-to-all pattern it increases linearly with the number of micro- 
modes. This quite surprising behavior not only contradicts our naive intuition 
that quantum correlations are averaged out in macroscopic systems as it em- 
phasizes that, the complexity of physical interactions between the microscopic 
constituents, can enhance the degree of macroscopic bipartite entanglement. 

However, there is a minimum interaction time needed before the two wave- 
packets start to be entangled, the birth time of entanglement, which increases 
with increasing temperature. 

Moreover, for the same interaction time and the same number of micro- 
modes, the degree of entanglement decreases with increasing temperature, be- 
coming zero at the critical temperature, where its first order derivative has a 
discontinuity. To associate this discontinuity with a possible quantum phase 
transition [27l [28] is not clear within our model. It would be interesting, in a 
further work, to use the techniques developed in solid state physics to relate the 
polarization and the temperature of the nonlinear medium where the waves in- 
teract, with their degree of entanglement, in order to find signatures of quantum 
phase transitions. 

Since this macroscopic entanglement is higher than the corresponding micro- 
scopic one and since the measuring apparatus only needs to detect the collective 
behavior of the wave-packets, it could be used as a resource to perform com- 
munication and computation tasks, at room temperature, in an efficient and 
technically more accessible way. 
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6 Appendix A: Microscopic evolution 



In this Appendix we solve the Heisenberg equations of motion for the Bose 
operators of the micromodes of the wave-packets. 

The solution of the Heisenberg equations of motion for di.fc and d\ _ fc in the 
pairwise interaction, was derived in |21j . and is given by 

Oi,k(*) = e _ ' i " 1 - fci (coshwt a w ,k + sinhiut a\ Q _ fe ) (22) 

a\_ k {t) = e tU2 - kt (smhwt a 10 , k + coshwt a|o,-fe) ( 23 ) 

The zero stands for initial time. 

For clarity of exposition let us rename the bosonic operators ai t k in the signal 
wave-packet 1 and 0,2,/ in the idler wave-packet 2 for the one-to-all interaction. 
The correspondence is the following: for modes of 1: a\— m — > Si;ai_ m +i — > 
a 2 ; a M „ -t a 2m +i and for modes of 2: a 2 ,- m -> b 2m+ 2; ■■■ a 2 , m hm+2- 
With this new notation the Hamiltonian H\ becomes 

2m+l 4m+2 2m+l 4m+2 

fli = Tw k ala k + + iw% J2 {e~' KQot+ip) albj-h.c.) 

fc=l /=2m+2 fe=l (=2m+2 

(24) 

The Heisenberg equations of motion for the vector operator Y = 
{&i,&2, —,Oam+ij4m+2'^2m+3) - ^L+2} T i under the Hamiltonian #l are 

= - lN[t )Y{t) (25) 

where iV(i) is 2n x 2n time dependent matrix. In order to obtain an autonomous 
system of equations we define new annhilation operators Z = F(t)Y where F(t) 
is a 2n x 2n diagonal matrix whose elements are 

f e iQ!t if i = j = l,...,2m + l 
Fij(t)=< er 1 ^ 1 if i = j = 2m + 2,..., 4m + 2 

The new bose operators obey the following autonomous system of linear differ- 
ential equations 

dZ(t) 



dt 

where 



MZ(t) (26) 



is the 2n x 2n matrix. A is a diagonal n x n matrix, whose elements are 



iS[m — (i — 1)] if i = j = 1, 2m + 1 
1 + 3 
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and B is the n x n constant matrix whose elements are 

Bij = w; = 1, 2m + 1) 



The solution of ( 26 1 is 

Z{t) = S*e- Dt S'- 1 Z(0) = G(t)Z(0) (28) 

Where D is the a diagonal matrix whose elements are the eigenvalues of the 
matrix M, and S is the matrix whose columns arc the eigenvectors of matrix 
M. Going back to the original bose operators Y(t), we obtain 

Y(t) = F _1 (t)G(t)y(0) = R(t)Y(0) (29) 

More explicitly, the time dependent operators cifc(i) and bj(t) in the one-to-all 
interaction are given by, 

2m+l 4m+2 

&k(t)= E r kj (t)& j0 + J2 r *i(*)*Jo ( 3 °) 

j=l 3=2m+l 
2m+l 4m+2 

&«(*) = E r «(*)a i0 + E ry(i)6} (31) 

j=l j'=2m+2 

in terms of the elements ry(i) of matrix i?(£). The subscript zero refers to initial 
time. 

The canonical commutation relations obeyed by Y imply that, 

2m+l 4m+2 2m+l 4m+2 

£hw(*)l a - E k«(*)l a = i ; EM*)I 2 - E hi(*)l 2 = -i 

j=l j=2m+2 j=1 j=2m+2 

(32) 

fc = 1, ...,2m + 1 and Z = 2m + 2, ...,4m + 2. 

7 Appendix B: Covariance matrix for the col- 
lective operators 

In this Appendix we compute the elements of the CM £ of the collective oper- 
ators. 

The collective time-dependent amplitudes in the one-to-all interaction are 

_^ 2m+l 2m+l 4m+2 

Mt)= /^==j E ^(t)= E "*j(*)^io+ E n i(*)*Jo ( 33 ) 

V^m+1 fc=1 j=2m+2 
4m+2 2m+l 4m+2 

£ f ffl = ^— rf E S K*)= E *i (*)&*)+ E %(*) s l-o (34) 

V^m+l ;=2m+2 J= 2m+2 
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where the coefficients 



^ 2m+l 2m+l 

= E r w (t); »,(*) = ^_ £ r w (t) (35) 

fe=l v k=l 

^ 4m+2 j 4m+2 

wo - ^g^r E ;.,('); «,(*) = 75^ E •*«) <»> 

/— 2m+2 2m+2 

are expressed in terms of the entries of the matrix R{t). The expressions for 
the elements of the symmetric collective covariant matrix for the initial thermal 
state, are easily computed 

£n(t) = E 22 (i) = + ! (37) 

2m+l 4m+2 

(Mt))= £(%) + l)K(t)| 2 + 2 %>l»i(*)l 2 (38) 

j=l ]=2m+2 

S 33 (*) = S44W = (N 2 (t)) + \ (39) 

2m+l 4m+2 

<#»(*)> = E %o|^)| 2 + E ("io + l)|«i(*)|2 (40) 

j=l ]=2m+2 

El 2 (t) = £34(t) = ( 41 ) 

{2m+l . 4m+2 . ~j 

Y J {n jQ + -)m j {t)t*{t)+ E (%+2K(iK( f )) ( 42 ) 
j=l " j=2m+2 " J 

S 24 (i) = -El 3 (t) (43) 



{2m+l ^ 4m+: 

E ("jo + 2 K (*)*,■(*)+ E 
j = l ~ j=1m- 



4m+2 _^ 

feo + ^Kf^iW)} ( 44 ) 

=2m+2 



= £i 4 (t) (45) 
The coefficients of the CM for the pairwise interaction are 

Sn(t) = E 2 2(t) = <JVi(t)> + ^ (46) 

(ATi(t)) = iVio cosh 2 wt + (iV 20 + 1) sinh 2 wt (47) 
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S 33 (i) - E44W - (N 2 (t)) + \ 



(N 2 (t)) = N w cosh 2 'wt + (N 20 + 1) sinh 2 'wt 
S 12 (t) = S 34 (i) = 



Sis(t) = i sinh(2u;t) cos(2u> i) ( ^10 + ^20 + 1) 



S 24 (i) = -E 13 (t) 



Su(t) = ^ sinh(2wi) sm(2cD t)(N w + N 20 + 1) 



£ 23 (t) = S i4(0 



Where 

Mo = <JV 1 (0))= X! 2^Tl ' ^0 = (^ 2 (0))= X! 



2m+l _ 4m+2 



+ 1 ' " ^ 2m + 1 

j=l j=2m+2 

are the initial mean number of photons in the wave-packets. 
For both interaction patterns 

det7(t) = -[E; 3 (t) + E? 4 (t)]<0 

For the initial vacuum state, n j0 = for any j, therefore N w = ft 
and dct j(t) = -sinh 2 wt < for t > 0. 
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G 2 



Figure 1: Graphs G\ and G2 represent the interaction between the microscopic 
modes of wave-packets 1 and 2 for Hamiltonians H\, panel (a), and H 2 , panel 
(b). Each vertex represents a microscopic mode of frequency ujij and each edge 
represents the interaction between mode (1, k) of wave 1 with mode (2, 1) of wave 
2. The coupling strength w is the same for all edges. 
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3.0 r 




Figure 2: Dimensionless BTE, Tg, over the temperature, for the dimensionless 
frequencies uj\ = uii/g = 200, W2 = ^2/9 = 400 and for different values of the 
dimensionless mismatch parameter y: y — (triangles), y — 0.5 (circles) and 
y = 0.9 (stars). 
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-r 



Figure 3: Dimensionless BTE, Tg in the case of the one-to-all interaction is 
displayed as a function of the temperature T for n modes in each wave-packet: 
n = 3 (crosses), n — 5 (inverted triangles) , n — 7 (stars), n — 9 (diamonds) 
and n — 11 (triangles). We assume the same dimensionless frequencies as in 
Fig. 2. 
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Figure 4: Dimensionless birth time of entanglement, T£ is plotted as a function 
of the number n of modes in each wave-packet for equilibrium temperatures: 
T = 30°K (stars), T = 300° if (triangles), and T = 3000° X (circles). We 
assume the same dimensionless frequencies as in Fig. 2. 




Figure 5: Logarithmic negativity £jj computed at the interaction dimensionless 
time Ti = 0.3324 is displayed as a function of the number n of micromodes 
for the initial vacuum state (triangles) and for a thermal state at equilibrium 
temperature T — 30°^ (circles) and for the one-to-all interaction. The same di- 
mensionless frequencies as in Fig. 2 are assumed. The solid lines that fit the log- 
arithmic negativity, are parallel straight lines whose equations are £tf — 0.963 n 
and Ejs/ — 0.963 (n — 5). 
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Figure 6: Logarithmic negativity £_\f, computed at the interaction dimensionless 
time Ti = 0.3324, is displayed as a function of the temperature T . Panel (a) n = 
1 and panel (b) n = 5, for the one-to-all interaction. The same dimensionless 
frequencies as in Fig. 2 are assumed. The critical temperatures are T c = 4.38°X 
for n = 1 and T c = 1200° if for n = 5. 
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